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Abstract 

Using the concept of cracking, we have explored the influence of density 
fluctuations on isotropic and anisotropic charged matter configurations in 
General Relativity with “barotropic” equations of state, P = P(p) and 
P± — P±{p) and a mass-charge relation Q = Q(p). We have refined the 
idea that density fluctuations affect physical variables and their gradients, 
i.e. the radial pressure and charge density gradients. It is found that not 
only anisotropic charged models could present cracking (or overturning), 
but also isotropic charged matter configurations could be affected by den¬ 
sity fluctuations. 


1 Introduction 

Despite that there exist some consensus about the charge neutrality of astro- 
physical objects, the debate dedicated to the influence of the electric charge 
on the structure and evolution of relativistic self-gravitating systems starts 
almost with the theory of General Relativity. Following the early contribu¬ 
tions of Rosseland and Eddington [1,2] there is a long list of works concerning 
charge and relativistic compact objects (almost twenty thousands entries in 
scholar.google.com and more than six thousands in last five years). There has 
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been a renewed interest in this subject emerging from the very ingenious theo¬ 
retical proposals for mechanisms allowing for the appearance of electric charge 
in strange stars with quark matter [3-9] and the consequences on the compact 
object structure of particular charge distributions regarding concepts such as 
global and local neutrality (see [10-13]). Obviously, all the above models are 
worthless if they are unstable against fluctuations of their physical variables and 
efforts to study of stability of neutral and charged matter configurations are as 
old as the modelling itself. 

One of the most common strategies to study the stability of selfgravitating 
fluids is to follow the evolution in time of the perturbations of the physical 
variables. In General Relativity this methology was developed in 1964, for 
pascalian fluids, by S. Chandrasekhar [14-16] and later it was extended to non- 
pascalian -i.e. anisotropic fluids having unequal radial and tangential pressures: 
P r ^ P±- by Hillebrandt and Steinmetz in 1976 [17] and, more recently, by Dev 
and Gleiser [18]. Despite this effort, Herrera and coworkers [19-21], following 
a qualitatively different point of view, introduce the concept of cracking (or 
overturning) to describe the behaviour of fluid distributions just after its depar¬ 
ture from equilibrium, when total non-vanishing radial forces of different signs 
appear within the system. They illustrated, under this approach, that some 
other possible additional pathologies of the fluid could be present in a matter 
distribution. Recently, [22,23], considering models described by “barotropic” 
equations of state P = P(p) and P± = P±(p), extends this criterion to isotropic 
and anisotropic spheres under local perturbations of the density that affect not 
only the physical variables but also their gradient. This extension leads to an 
straight forward generalisation of the cracking criterium for anisotropic charged 
fluids which is developed in the present paper where we depart from the con¬ 
jecture that the charge and the mass of a fluid are related by a “barotropic 
equation of state”, Q = Q(p). Following the same scheme developed in [22,23] 
apply this idea to non constant (local) density fluctuations that can affect the 
gradient of pressure and charge density considering the local fluctuations, i.e. 
Sp = Sp(r), represented by any function of compact support defined in a closed 
interval A r <C 1. 

This work is organised as follows. Section 2 will present the context & 
notation. In Section 3 we describe the concept of cracking for selfgravitating 
charged matter configurations presenting the effects of local density fluctuations 
on the force distributions within a matter distribution. Section 4 illustrates these 
effects by some examples. Finally, some conclusions are displayed in Section 5. 


2 Hydrostatic equilibrium of charged spheres 

Let us consider an static spherically symmetric metric, described by the line- 
element 

ds 2 = e 2 " (r) d t 2 - e 2A(r) dr 2 - r 2 (d<9 2 + sin 2 6 dc/) 2 ) , (1) 
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with a matter component of the energy-momentum tensor that corresponds to 
an anisotropic fluid 

= (,p + P±)u IJ ,u l ,-P±g liU + (P-P±)v IJt v v +- ^g a pF a pF a0 ^J , 

( 2 ) 

where p is the energy density, P the radial pressure, P± the tangential pressure, 

= (e",0,0,0), and = (0, -e A , 0,0) , (3) 

are respectively the unitary velocity vector of a co-movil observer and a unitary 
normal vector to it. The electromagnetic field tensor F /t „ can be written in 
terms of the electromagnetic potential A 

Fin* = i ( 4 ) 


therefore, the Maxwell electromagnetic field equations are 

{^~gF^) v = An^r, (5) 


with = an, ^ the current four vector, and a the charge density. 

Since we want to consider an electric charged sphere, the potential vector is 
of the form 

Ap = (A), 0, 0,0) , (6) 

then the only components of the electromagnetic tensor are For and — Fio, and 
the Maxwell field equations (5) can be written like 

- (e- (l/+A) r 2 F 0 i) = 47rr 2 e A cr . (7) 

Defining charge Q(r) as 


Q{r) = f 47rr 2 (Te A dr , 

Jo 


we can express the components of the electromagnetic tensor Foi as 

F _ Q( r ) -\+V 

foi —-5—e 


and the non-null components of the energy-momentum tensor, F^„, as 


( 8 ) 

(9) 


T% = P+ ^e- ( - A) , (10) 

87 r 

r 1 ! = -p + £k e - (y+A) , ( 11 ) 

87T 

T 2 2 = T 3 3 = -P ± - ^ke“ (l/+A) . (12) 
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Following [24], we calculate r —component of the conservation equation, 


1 ;„ = 0 


( 13 ) 


and obtain 

dP 1 Q dQ | (P + P ) ( 4 ^ rP + ~ ft) 2 (P ± - P) _ Q 

dr 4 7rr 4 dr 0- ~ ~— + ^-) r 

which is the generalization of the Tolman-Oppenheimer-Volkoff equation of hy¬ 
drostatic equilibrium for a charged matter distribution. Since we are consid¬ 
ering a static charged sphere, the exterior region is described by the Reissner- 
Norsdstrom metric 

ds 2 = (V — + -SI) d* 2 - T - r 2 (dd 2 + sin 2 9 d cf) , (15) 

\ r r 2 J i_ 2 M_i_Qo v ' v ' 

r ' r 2 

where Q o is the total charge and M is the total mass of the sphere. The in¬ 
ner region is matched to the Reissner-Norsdstrom metric through a null radial 
pressure hypersurface of radius tq. 


3 Cracking of anisotropic charged spheres 


In this section, following our previous work [22,23], we extend the concept of 
cracking for isotropic and anisotropic charged matter configurations. 

As we have stressed before this concept was introduced to describe the be¬ 
haviour of fluid distributions just after they experiment a perturbation in any of 
its state variables and departs from equilibrium. Thus, the r—component of the 
conservation equation (14) which represents radial forces may be different than 
zero and radial forces may appear on the system which could be represented by 

d P 1 Q dQ (p + P)(47rrP+^-^) 2 (P ± - P) 

dr 4 7rr 4 dr p _ Zm _|_ Q^.'j r 

Herrera and collaborator [19-21] had shown that there exist a relation between 
the force distribution and the expansion written as 


e x (p + P r ) 

e u/2 r 2 


dr e y/2 r 2 


d0 

ds 


(17) 


where 0 represents the expansion. Therefore if there is a change of sign on 
the force for some radial coordinate within the configuration, it affect also the 
expansion and we say that cracking is produced. 

Now let us consider an anisotropic charge sphere in hydrostatic equilibrium, 
described by equations of state of the form P = P(p), P± = P± (p ), and that 
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hold a charge-mass relation of the form Q = Q(p). Then, we will consider a 
local density perturbation 

5p = Sp(r) , (18) 

described by any function of compact support Sp = Sp{r) defined in a closed 
interval A << 1, which generate fluctuations in mass, radial and tangential 
pressures, and pressure gradient that can be represented as 

dP 

P(p + Sp) ft P(p) + SP = P(p) + — Sp = P(p) + v 2 Sp , (19) 


Pj_(p + Sp) ft Pj_(p) + SP± = P±(p) + ——Sp = P±(p) + v 2 ± Sp , (20) 

d p 


P'(p + 5p) 


P'(p) + SP' = P'(p) + ^Sp 
dp 


= P’(P) + 


/ 2\/ , 2 P 
{v) + V — 


Sp , 


( 21 ) 


Q(p + Sp) ft Q(p) + SQ = Q(p) + ~Y~Sp = Q(p) + %-Sp , (22) 

d P p' 

dO' Q" 

Q'{p + Sp) « Q\p) + Sp = Q'(p) + \8 P , (23) 

op p' 


m(p + Sp) ft m(p) + 8{m) 


, , dm 4irr 2 p 

m(p) + 1 Sp = m(p) H- —Sp . 

op p' 


(24) 


Where primes denote radial differentiation while 


v 
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dP 2 _ d P ± 
dp ’ v± dp 


m = 47t / p(r)r 2 dr , 

J o 


(25) 


are the radial sound speed, tangential sound speed and mass, respectively. 

Now, we can, formally, expand equation (16), and obtain the radial force 
that appears due to the perturbation of density 


n^n a (p,m,pP , ,P±,Q,Q') + Sn, (26) 


with 


&R- dn : ondn sn , dn dn^ on 
67Z -~d^ 6p+ d^ Sm+ dP SP+ d^ 6P+ WI SP± + dQ 6Q+ W Q ’ ( } 


and 


dP dP' 

Po(p, m, P, P 1 , Pj_, Q, Q') = 0 , 


(28) 
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since the system is initially in equilibrium. By using (19) it can be shown that 


<57? = Sp 


dn dn 

dp + dm 


47rr 2 p 
P' 


dn 


dn 


-v + 
dP dP' 


dn 


dn 


dP± V± + dQ 


yy +v 7 ) 

<n .wfoax 

P' ) dQ' \p' )_ 


where 

dn 4irr 4 P + mr — Q 2 dn {p + P)(r 2 — Q 2 + 8nr 4 P) 

dp r(r 2 - 2mr + Q 2 ) ’ dm ( r 2 - 2mr + Q 2 ) 2 

97?. 1 dQ 2 (p + -P) (47rr 4 P + mr - Q 2 + l) <2 

dQ 4 ttt 4 dr r (r 2 — 2mr + Q 2 ) 


(29) 


(30) 


(31) 


dn 8nr 4 P — 3mr + Q 2 + 47t r 4 p + 2r 2 

dP r(r 2 — 2 mr + Q 2 ) 


(32) 


dn _ 1 dn _ 2 

dQ 1 4wr 4 ’ dP± r 



(33) 


As in [22,23], if <57? does not change its sign we will consider that the con¬ 
figuration is potentially stable [21,25]. Since this time there are terms of the 
charge, and the charge gradient, is not simple to identify which terms are the 
ones that may contribute to the change of sign of the radial force; then the whole 
expression must be taken into consideration. In the next section we will analyse 
this term in some models in order to identify instabilities under this criterion, 
studying how the terms Ri determine the outcome of the perturbation 

~ SR ~~~~~~~ 

R — — = -Ri + R2 + R3 + R4 + R5 R7 + Rs j ( 34 ) 

Sp 



dn 

dn 47rr 2 o 

- dn n 

- dn , 

(35) 

Ri = 

~dp’ 

H2 dm 

P' 

^3 = , 

R ‘-ap ± v ^ 

i?4 = 

dn 

ap 

r 0 " 1 

(V ) +V — 

[ P \ 

5 

p dn Q' 

Ke — > 
dQ p' 

p dn Q" 

7 dQ> p' ' 

(36) 


4 Analysis of charged models 

In this section we illustrate, through examples, the three possible cases -cracking, 
overturning and no cracking at all- for isotropic and anisotropic fluid relativistic 
spheres. 
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4.1 Model 1: Isotropic charged sphere 

We analyse an isotropic charged sphere model proposed in [26], equations (3.17) 
and (3.18), described by the density, pressure, and charge distribution: 


P = P r o 


5 f 2 — 3 — 2>a±q 2 + 3/xa 4 — 3 a 2 q 2 q 2 

87 tq:| 47 T r 2 ’ 


(37) 


P = Prl = 


+ 


2 yja\ + (1 + a±q 2 - pot\ + a\q 2 )f 2 - r 4 - a 2 q 2 
8na 2 

l + «4 q 2 {q 2 — n)a\-r 2 q 2 

87 ra| + 8irr 2 ’ 


(38) 


Q = Qr%=qr, 


(39) 


P± = P, 

with 

1 (12/x — 12 — q 4 )(fJ? — 4g 2 /r + 4g 4 ) 

“ 4 ~ ~ 3 (// - 2<? 2 )Mm - 2 9 2 ) 2 (/r 2 - 4g 2 At + 4g 4 ) 2 ) 1 /3 

3 (p - 2g 2 ) ’ 


(40) 


(41) 


x (a 2 ^-2q 2 ) 2 (^-4q 2 „ + ^) 2 y /3 
3 (/z — 2q 2 ) {fj? — 4q 2 [i + 4 q 4 ) 


(42) 


a 2 = 18 q 2 qL + 36g 2 — q 6 — 108g 4 + 6\/3ai , (43) 

ai = [16/x 3 (—48 - q 4 )fi 2 + (48 + 20<? 4 - 3 6q 6 )fi - 16 

+2g 1 0 + 8 q 4 + 107 q 8 - 72g 6 ] 1/2 . (44) 

Now, substituting the expressions for the density, charge and pressures, using 
the parameters /i = 0.1 and q = 0.1. In figure 1 we plot i?, of equation (34), and 
the terms Ri that compose it separately, equations (35) to (36). This model 
exhibits cracking since R changes sign around r = 0.7, this indicates that around 
that region the model could be potentially unstable under perturbations of the 
density. 
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Figure 1: On the left, the force distribution, 1Z = 51Z/8p for the isotropic 
charged sphere proposed in [26], using the parameters p = 0.1 and q = 0.1. 
Since the force changes sign around f = 0.7, we say that the configuration could 
be potentially unstable under the perturbation of density 8p the region near that 
coordinate. On the right plate, all the terms that compose the force distribution 
IZi are sketched . 


4.2 Model 2: Anisotropic charged sphere I. 

Now, we analyse an anisotropic charged model proposed in [27], described by 


2 1 

P = P r o = o 


n + f 


8 7T (1 + 2n — C) r 2 


(45) 


P = Pr z n = - 


1 + f 2 


8 7r (1 + 2n — C) r 2 


(46) 


P± = P±r 2 0 = CP , 


(47) 


Q = Qr 0 


r 2 (1 + n — 2C — r 2 (l — C)) 


2 + 4n - 20 


(48) 


This model represents a matter distribution, admitting a one parameter group 
of conformal motions with a linear equation of state. 

Now, substituting the expressions for the density, charge and pressures, using 
the parameters O = (n— 1)/2 — 0.1 and n = 2. In figure 2 we plot R, of equation 
(34), and the terms Ri that compose it separately, equations (35) to (36). Now, 
as well as the previous example, this model is potentially unstable since R 
changes its sign around f = 0.6. 
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Figure 2: On the left, the force distribution, 1Z = 51Z/5p for an anisotropic 
charged sphere proposed in [27], using the parameters C = (n — l)/2 — 0.1 and 
n = 2. Since the force changes sign around f = 0.6, we say that the configuration 
could be potentially unstable under the perturbation of density Sp around the 
region near that radial coordinate. Again, on the right plate all the terms that 
compose the force distribution IZi. 


4.3 Model 3: Anisotropic charged sphere II 

Now, let us consider an anisotropic charged model proposed also in [28], equa¬ 
tions (125), (126) and (127), described by the dimensionless equations 


p = prl = 


(3Ci - 1) r (1-Ci) 
87Tf 2 [f 2 (3Ci - 1) 


+ 1 + 




(49) 




2 ) (1 - a 2 r 2 ) 


(50) 


P± 


P±r 2 


(1 ~ 2 CQ 

87rr 2 


+ P, 


(51) 


Q 


Q \/3aCi — af 3 
To 2 


(52) 


Now, substituting the expressions for the density, charge and pressures, using 
the parameters C\ = 0.4 and a = 1, in figure 3 we plot R , of equation (34), 
and the terms Ri that compose it separately, equations (35) to (36). As well as 
the previous example, this model has cracking instability since 1Z changes sign 
around r Ri 0.997. 
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Figure 3: On the left plate, the force distribution, 1Z = STZ/Sp for an anisotropic 
charged sphere [28], using the parameters C\ = 0.4 and a = 1 . Since the force 
does not change sign, no cracking occurs and the configuration is stable under 
the perturbation of density Sp. On the right, the terms that compose the force 
distribution IZi. 

5 Remarks and conclusions 

We have shown that isotropic and anisotropic charged matter configurations 
can present cracking and overturning when non constant (local) density fluctu¬ 
ations are considered and affect the state variables and their gradient. As we 
have pointed out, the density perturbations we assume are local, Sp = Sp(r ) 
-represented by any function of compact support defined in a closed interval 
Ar -C 1- within a fluid configuration governed by barothropic equations of 
state: P = P(p), P± = P±(p) and a mass-charge relation Q = Q{p). 

The idea of cracking was originally conceived by Herrera [19] to describe 
the behaviour of fluid distributions just after their departure from equilibrium: 
fluid elements, at both sides of the cracking point are accelerated with respect 
to each other by -independent and simultaneous- perturbations in energy den¬ 
sity and anisotropy [20,21], In their approach they considered independent and 
simultaneous perturbations in density and anisotropy that lead cracking (or 
overturning) within the configuration. Later, Abreu, Hernandez and Nunez [25] 
shown how constant (global) density perturbations could generate cracking on 
anisotropic relativistic fluids. In this study, constant density fluctuations affect 
mass, radial and tangential pressure, but leave unperturbed the pressure gra¬ 
dient; again only anisotropic distributions can exhibit cracking or overturning. 
This idea was extended to the charged case in [29] where some preliminaries 
results were presented. In a recent work [22,23], we extend this criterium -of 
density perturbation affecting the state variables and their gradient- and showed 
that cracking (or overturning) could also be present for isotropic and anisotropic 
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spheres. This extension leads to an straight forward generalisation of the crack¬ 
ing criterium for isotropic/anisotropic charged fluids developed in the present 
paper. 

It is worth to be mentioned that, the concept of cracking illustrate that 
some pathology could be present for a particular equation of state and refers 
only to the tendency of the configuration to split (or to compress) at a par¬ 
ticular point within the distribution but not to collapse or to expand. The 
cracking, overturning, expansion or collapse, has to be established from the in¬ 
tegration of the full set of Einstein equations, but it should be clear that the 
occurrence of these phenomena could drastically alter the subsequent evolution 
of the system. If within a particular configuration no cracking (or overturn¬ 
ing) is present, we could identify it as potentially stable (not absolutely stable), 
because other types of perturbations could lead to its expansions or collapse. 
Within a relativistic matter configuration, in principle, it is not clear which of 
the two scenarios is more likely to occur: the simultaneous two-perturbation 
original scenario of Herrera and coworkers or the only density-driven framework 
but surely both generate instabilities on relativistic matter configurations that 
should be evaluated. 
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